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Dewey  H.  Hodges. 3  The  subject  paper  deals  with  mathematical 
modeling  of  the  nonlinear  behavior  of  beams.  A  set  of  equations  is 
derived  and  used  to  investigate  the  static  behavior  of  a  slender  can¬ 
tilever  beam  loaded  transversely  at  the  free  end.  Since  the  loading  is 
not  necessarily  along  the  principal  axes,  the  principal  bending  de¬ 
flections  and  torsion  are  coupled  elastically.  In  a  strictly  linear  theory 
there  would  be  no  torsion,  thus  the  presence  of  torsion  is,  itself,  a 
manifestation  of  nonlinear  behavior.  The  agreement  obtained  with 
experimental  data  achieved  in  the  subject  paper  is  excellent  and 
confirms  that  one  may  apply  the  geometric  nonlinear  theory  of  elas¬ 
ticity  to  beams  with  confidence  provided  an  adequate  degree  of 
nonlinearity  is  retained  in  the  mathematical  model  to  account  for  large 
deflections. 

The  same  problem  was  treated  by  Dowell,  et  a!.,  in  |1 1.  The  ana¬ 
lytical  results  of  1 1 1  and  those  of  the  subject  paper  agree  quite  well  as 
long  as  ttt.  the  bending  deflection  in  the  plane  of  greatest  flexural 
flexibility,  remains  relatively  small  compared  to  the  beam  length.  In 
fact,  analytical  results  of  1 1 1  depart  from  those  of  the  subject  paper 
only  when  w  ceases  to  be  a  small  fraction  of  the  beam  length.  An  ex¬ 
amination  of  the  w  equation  (.13 b)  in  the  subject  paper  reveals  that 
in  addition  to  linear  and  second -degree  terms,  several  terms  of  third 
degree  in  the  deflections  appear.  These  terms  are  retained  in  the  au¬ 
thors'  original  derivation  of  these  equations  [2|  for  the  special  case 
when  hiJlsi  is  large  compared  to  unity.  Although  not  stated  in  the 
subject  paper,  these  third-degree  terms  constitute  the  only  difference 
between  the  equations  of  |1|  and  the  subject  paper.  Thus  the  im¬ 
proved  agreement  reported  in  the  subject  paper  must  be  due  to  the 
presence  of  these  third-degree  terms  and  not  to  the  reasons  indicated 
in  the  paper.  The  two  main  reasons  given  for  the  improved  agreement 
shown  in  the  subject  paper  are  now  examined. 

The  first  reason  is  mentioned  on  Page  162.  Column  1 ,  where  the 
authors  attribute  the  accuracy  of  their  results  to  "a  more  careful  and 
consistent"  derivation  than  in  |1|.  The  original  derivation  of  the  au¬ 
thors’  equations  [2|  concludes  with  a  set  of  nonlinear  equations  of 
second  degree  in  bending  and  torsion  deflections.  These  equations 
are  intended  for  applications  to  rotating  blades  and  are  based  on  the 
assumption  that  bending  and  torsion  rigidities  are  of  the  same  order 
of  magnitude.  For  the  special  case  where  /r.//: n  and  are  large 

compared  to  unity,  the  authors  endeavored  to  retain  third  degree 
terms  whenever  they  are  multiplied  by  a  large  coefficient  in  the 
equations.  It  is  not  clear  that  appropriate  measures  were  taken. 


1  By  A.  Rosen  and  P  Friedmann,  and  published  in  the  March,  1979,  issue 
of  the  ASME  Journal  op  Appi.if.d  Mechanics,  Vol.  46,  pp  l6i  168. 
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however,  to  ensure  that  all  third -degree  terms  multiplying  the  large 
coefficients  were  retained.  In  fact,  only  terms  through  second  degree 
were  retained  in  the  bending  curvature  expressions,  and  it  was  never 
demonstrated  that  third-degree  terms  from  these  expressions  would 
not  appear  in  the  final  equations.  Moreover,  the  final  equations  in  the 
subject  paper,  because  of  the  particular  third-degree  terms  retained, 
do  not  have  a  self-adjoint  structural  operator.  While  it  is  acknowl¬ 
edged  that  the  authors  deserve  credit  for  adding  appropriate 
higher-degree  terms  to  the  equations  to  improve  the  correlation  with 
experimental  data,  this,  in  itself,  does  not  necessarily  imply  a  higher 
degree  of  care  or  consistency. 

The  second  reason  given  by  the  authors  appears  in  the  C oncluding 
Remarks  section:  “The  superior  agreement  obtained  with  the  present 
equations  is  due  to  differences  between  the  final  equations  of  equi¬ 
librium  used  here  and  those  given  in  |3],  These  differences  have  been 
discussed  with  considerable  detail  in  [2]."  (Reference  numbers  refer 
to  those  of  this  discussion).  The  first  sentence,  by  itself,  appears  to 
be  referring  to  the  third -degree  terms  retained  in  equation  (33b  >  be¬ 
cause  the  ratio  Iw/ha  is  large  compared  to  unity  for  the  particular 
beam  considered.  This  is  not  the  intended  meaning  according  to  the 
second  sentence,  however,  because  differences  related  to  these  third 
degree  terms  are  not  discussed  in  |2|.  The  differences  between  the 
equations  of  [2, 3j  that  are  discussed  in  |2)  relate  to  the  second -degree 
equations  only  and  have  been  recently  clarified  in  |4,  5).  The  final 
second-degree  nonlinear  equations  of  |2,  3),  while  differing  slightly 
for  pretwisted  beams  |4],  are  actually  equivalent  for  the  present  case 
of  a  nonrotating  beam  without  pretwist  |5).  Therefore,  the  improved 
agreement  reported  in  the  subject  paper  does  not  prove  that  the 
equations  for  rotating  blade  applications  derived  in  |2)  are  more 
"reliable"  nor  that  they  "can  be  used  with  confidence"  any  more  than 
those  of  |3|.  The  results  simply  prove  that  the  equations  of  the  subject 
paper  are  more  accurate  than  those  of  1 1 1  when  applied  to  problems 
in  which  the  bending  deflection  ir  may  exceed  a  small  fraction  of  the 
blade  length  and  the  ratio  Ijx/lm  is  large  compared  to  unity.  These 
cases  are  not  treated  in  |l )  in  which  the  squares  of  bending  slopes  were 
assumed  to  be  negligible  with  respect  to  unity  and  the  equal  ions  were 
taken  directly  from  |3]  where  1?>  and  / were  assumed  to  Ire  the  same 
order  of  magnitude  |1,  pp.  8. 9|. 
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